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We discuss the non-linear gravitational interactions in the Randall-Sundrum single brane model. 
If we naively write down the 4-dimensional effective action integrating over the fifth dimension 
with the aid of the decomposition with respect to eigen modes of 4-dimensional d'Alembertian, the 
Kaluza-Klein mode coupling seems to be ill-defined. We carefully analyze second order perturbations 
of the gravitational field induced on the 3-brane under the assumption of the static and axial- 
symmetric 5-dimensional metric. It is shown that there remains no pathological feature in the 
Kaluza-Klein mode coupling after the summation over all different mass modes. Furthermore, 
the leading Kaluza-Klein corrections are shown to be sufficiently suppressed in comparison with the 
leading order term which is obtained by the zero mode truncation. We confirm that the 4-dimensional 
Einstein gravity is approximately recovered on the 3-brane up to second order perturbations. 
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I. INTRODUCTION 



The possibility that our 4-dimensional universe is embedded as a brane in a higher dimensional spacetime has been 
extensively discussed recently as the brane world scenario pi— HB . In particular, Randall and Sundrum introduced 



' attractive models whose background bulk spacetime is 5-dimensional anti-de Sitter (AdS^) spacetime. These models 
. suggest the possibility of the existence of extra dimensions in the non-trivial form and the possible explanation of the 
^\ ' large hierarchy between the Planck scale and the electroweak scale || Ij . 

The behavior of the gravity in the Randall-Sundrum (RS) models has been investigated by many authors [p[dl4[ . 
• For the RS single brane model, which is the model with the positive tension brane alone, the extension of extra 
dimension is infinite. Nevertheless, the gravity on the 3-brane generated by the matter field, which is confined on the 

3- brane, approximately coincides with the 4-dimensional Einstein gravity 

For the RS two-brane model whose bulk space is bounded by the positive and the negative tension branes, the 

4- dimensional Einstein gravity can be recovered on both branes under the approximation of the zero mode truncation 
H when we take the stabilization mechanism of the distance between the two branes into consideration jl^ [l7| . 

Although a large number of studies have been made on linear perturbations of the metric |lj |23|] , little is known 
about the non-linear or non-perturbative feature of the gravity p4|-|30| . There is a worry about the RS single brane 
model. In this model, it seems that the non-linear gravitational interaction between Kaluza-Klein (KK) modes becomes 
strong and it diverges as we move far away from the brane. A more precise argument is that the 4-dimensional effective 
. ^ action including interaction terms diverges when we attempt to write it down by integrating over the dependence 
. on the 5th direction by using the decomposition of metric perturbations in terms of eigen modes of 4-dimensional 
d'Alembertian jjj. This fact indicates that we cannot construct a 4-dimensional effective action for this model in an 
5^ , usual sense. 

The above discussion is based on the analysis using the decomposition of the mass eigen modes, which is referred 
to as the mode-by-mode analysis. However, it was demonstrated that the mode-by-mode analysis is insufficient to 
deal with metric perturbations and it is necessary to take into account the contributions from all the KK modes when 
we discuss the regularity of linear metric perturbations at a point far from the 3-brane jnj. Thus the pathological 
behavior in the non- linear interaction can be expected to be an artifact solely due to the mode- by- mode analysis, 
although it has not been proved yet. Even if the pathological behavior might be fictitious, there is another question 
if the 4-dimensional Einstein gravity is recovered when we proceed to higher order perturbations? 

The purpose of this paper is to give a partial answer to the above two questions. To study the non-linear behavior 
of the gravity, we investigate second order perturbations in the context of the RS single brane model. To simplify the 
analysis, we consider the static and axisymmetric configuration, which means that the 3-brane metric is spherically 
symmetric. Following the method developed by Garriga and Tanaka Q, we confirm that there is no pathological 
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feature in the KK mode coupling if we sum up all the mass eigen modes, and that the contribution due to the KK 
mode coupling is sufficiently small compared with the leading order terms, which are obtained by the zero mode 
truncation. Furthermore, the result obtained by the zero mode truncation exactly agrees with the one predicted by 
the 4-dimensional Einstein gravity, n 

This paper is organized as follows. In the next section we first derive the second order perturbation equations in 



the 5-dimensional bulk, and discuss the gauge transformation as well as the boundary condition. In Sec. Ill, after 



giving brief summary on the behavior of the 5-dimcnsional Green function, which is need to solve the perturbation 
equations, we review the results for linear perturbations to give their explicit expressions in the notation of the present 
paper. In Sec. IV, we analyze second order perturbations of the metric induced on the 3-brane. First it is shown that 
the 4-dimensional Einstein gravity is recovered in the approximation of the zero mode truncation, and then we prove 
that the remaining contribution due to the Kaluza-Klein mode coupling can be neglected. Section [v| is devoted to 
summary. 



II. PERTURBATION EQUATIONS IN THE RS MODEL 

The brane world model proposed by Randall and Sundrum is composed of the 5-dimensional AdS space 

ds 2 = g ab dx a dx b = a 2 (y)j] tiv dx^dx v + dy 2 , (2.1) 

with a single positive tension (cr > 0) 3-brane located at y = 0. Here a(y) = e~\ y ^ e is the warp factor, and £ is the 
curvature radius of AdS$. We have denoted the 4-dimensional Minkowski metric as rj^. The cosmological constant 
on the bulk and the tension of the 3-brane are, respectively, related to the curvature radius £ as A = — 6£~ 2 and 
as a = 3/At:£Gc i , where G5 is 5-dimensional Newton's constant. The relation between the 5-dimensional and the 
4-dimensional Newton's constants is given by G5 = IG±. We also use the notation k$ = SttG^. Ordinary matter field 
is supposed to be localized on the brane. 

In this paper, we investigate second order perturbations of this model induced by non-relativistic matter on the 
brane whose energy momentum tensor is given by the perfect fluid form 

T; = dmg{- p.P.P.P). (2.2) 

To simplify the analysis, we restrict our consideration to the static and axisymmetric spacetime whose axis of symmetry 
lies along y-direction. Namely, the 3-brane metric is spherically symmetric. 

We denote the perturbed metric by g a b = gab + ftab- At the level of linear perturbations, it is advantageous to use 
the Randall-Sundrum (RS) gauge defined by 

h"yy — ft/xy — 0, (2.3) 

ft V = 0, ft% = 0, (2.4) 
because the linear perturbation equations in this gauge take the simple form B 



r 2 D (4) + d 2 - Ar 2 



ft M „ = 0, (2.5) 



and all components are decou pled . However, whe n w e consider second order perturbations, we cannot imp ose the 



transverse-traceless condition ( |2.4|) in addition to (2J5). Hence, we need to abandon to impose the condition (2.4) in 



second order perturbations. As a consequence, the second order perturbation equations are inevitably coupled. 



1 After we had submitted the previous version of this paper, we had the existence of Ref. [g^J pointed out. In this reference, 
the second order metric perturbation in the region far from the star was investigated using the truncation of the first order 
metric perturbation at the leading order for the expansion with respect to the distance from the star. Since the Einstein 
equation was not solved in the whole region of the spacetime in their treatment, there remains an ambiguity to add some 
metric perturbations which satisfy the homogeneous linear perturbation equation only in the asymptotic region. Nevertheless, 
one may be able to prove that the dominant part of this remaining ambiguity can be absorbed by the redefinition of the mass 
parameter, and actually our results show that this is the case. Hence, the recovery of the 4-d Einstein gravity in the region 
far from the star is to be credited to Ref. |m| . On the other hand, what we show in the present paper is the recovery of the 
4-d Einstein gravity throughout the whole spacetime including the inside of the star without assuming a specific radial matter 
distribution. 
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Here, instead of requiring (2.4), we start with the assumption of the diagonal form of metric 



ds 2 = a 2 



-e A ^dt 2 + e B ^dr 2 + e c{r ^r 2 (d6 2 + 



sin 



) +dy 2 



(2.6) 



which is manifestly compatible with the gauge condition (2.3), and does not loose the generality under the restriction 



to the stat ic a nd axial-symmetric case. Furthermore, we can expect that this assumption is also compatible wit h the 
condition (2.4) at the linear order according to the result obtained in Ref. S. For the assumed metric form (2.6), the 
conditions (2.4) at the linear level, respectively, become 



A™ + 2C« = 0. 



and 



d r {r 2 B {1) )-2rC {1) =0, 
where we have expanded A, B and C to the second order like 

A(r,y)= ]T A {J Hr,y). 



(2.7) 



(2. 



(2.9) 



.7=1,2 



Hereafter, we neglect higher order terms without mentioning it. 

Before we start to solve the 5-dimensional Einstein equation, we would like to mention the boundary condition 
at y — + oo. An important point which we must mention here is that we are to find a solution which is regular at 
y — > oo. If we allow the violation of regularity at the infinity, the dynamics of the RS brane world model is not 
uniquely determined. Then, such a model cannot be a candidate of the model that describes our universe. Hence, 
we require that the metric converges to the AdS^ asymptotically. To guarantee this asymptotic condition, we simply 
require that the metric functions A, B and C to go to at y — ► oo. With this choice, the residual gauge degrees of 
freedom are completely fixed. Under this constraint, it is still possible to extend the coordinates throughout the bulk 
maintaining the diagonal form of the metric. We extensively refer to this gauge choice as the RS gauge because we 
see later that the first order quantities A^',B^ and C^ 1 ' in this gauge satisfy the both cond itio ns (2.7) and ( |2.8| ) (, 
although the quantities at the second order do not satisfy the transverse-traceless condition (2.4)). 



A. 5D Einstein Equations in the bulk 



We consider the 5-dimensional Einstein equations and derive the perturbation equations up to the second order. 
Since the trace of metric perturbations vanishes in the RS gauge at the linear order, it would be convenient to introduce 



2^ (J) =^ (J) +/3 (J) +2C( J ). 
By using this quantity, the trace of metric perturbations is expressed as 



rv = 2 E ^ (J) 



.7=1,2 



(A^) 2 + (B^) 2 + 2(C^) 2 



(2.10) 



(2.11) 



Hence the traceless condition at the first order is simply given by = 0. 

The 5-dimensional vacuum Einstein equation with the cosmological term is equivalent to the following set of 
equations for the Ricci tensor, 



4 



,7=1,2 



,7=1,2 



(2.12) 
(2.13) 
(2.14) 



.7=1,2 



where A = X)i=i ®1 an d we have defined 
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S(r, y) = ~d r (t*A$>BM) + a 2 ^A% + , (2.15) 

Q^y) ee \ [A<$A$ +B«BW + 2CWcW) . (2.16) 

Note that no other equations for the remaining components are independent. 

First we consider the trace part i])^ J \ which can be evaluated by integrating ( 2.12 ). From the requirement of the 
boundary condition, ipW must go to at y — > oo. Thus, we obtain 



^ (J) (r,»)=-e M / dy"e 2y "^ 



dy e 



i p -ly'/ 



iyy 



(2.17) 



where we have introduced a symbol e' J ) that is defined by = and = 1 to r epre sent the first and the second 
order equations in a single expression. As is anticipated above, we can see from Eq. (2.17) that the traceless condition 
at the linear order is actually satisfied, while that at the second order cannot be imposed any longer in general. 
Let us now turn to B^ J K Integrating ( [2.14 ), we obtain 



£( J )(r,j,) = -l J r 2 A^dr + ^e^ (J drr 2 (r^), r 



+ dy dr r Q ry 



(2.18) 



As before, the integrat ion constant is fixed by the boundary condition at y — ► oo. Now, combined with Eq. ( p.7j ), it 
is easy to see that Eq.(2.8) holds at the linear order. Hence, it is confirmed that our choice of gauge is equivalent to 
the RS gauge. 

So far, we obtained the relations between the metric functions A^ J \ B^ andC^' 7 ^. Substituting these into Eq. (2.13), 
we obtain a single equation for A^ J \ 



£[fl 2 iW]=eW[2a 2 rV, 9 -S] 



(2.19) 



The source terms are absent at the linear order, as they arc composed of the quadratic in the linear order quantities. 
Since ?/>W = 0, S defined in Eq.( [2.15| ) simplifies as 



S(r,y) = ^d r (r 2 A^B^ 



(2.20) 



B. Boundary condition 



In the previous subsection, we have obtained the master equation for the metric functions in the bulk up to the 
second order. To solve this equation, we need to know the boundary condition to be imposed on the 3-brane. The 
boundary condition on the 3-brane is specified by the Israel's junction condition pi] . However, in the RS gauge 
defined above, the y-constant surface is, by construction, chosen so that the metric functions go to at y — > oo. After 
fixing the y-constant surface for large y, coordinates are extended to the region near the 3-brane. Therefore, we can 
no longer expect that the location of the 3-brane coincides with the y = surface in general ||. In such coordinates, 
the junction condition is not so trivial. Thus it will be convenient to introduce other coordinates x a , in which the 
location of the 3-brane stays at y — but the metric form is still kept diagonal. We associate a bar with the quantities 
written in these coordinates, say, A. We denote this choice of coordinates by the Gaussian normal (GN) gauge. By 
construction, h^ly—o gives the 4-metric induced on the 3-brane. 

The junction condition on the 3-brane is simply written in the GN gauge as ^ 



2 One may think that this equation should have been expressed as 



p»(x)dy(a-\y)h^(x)) = -K 5 (T\ l -^\T\yx). (at y = 0+), 

because this is a relation for the quantities in the GN gauge. However, f(x)\ y= o = is obviously equivalent to f(x)\y= o — 
for any function /. To stress this point, we always use the coordinates without bar hereafter. The only exception is Eq. (2.24). 
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^(x)d y (a- 2 (y)h^(x)) = -«5 ( T\ - U\ T\ I (x). (at y = 0+) 



(2.21) 



As mentioned earlier, we assume the energy momentum tensor of the perfect fluid form (2.2). Then, the 4-dimcnsional 
energy-momentum conservation T v M; „ — becomes 



(p(r) + P(r))d r .4 (1) (r, 0) + 2 d r P(r) = , 



(2.22) 



and hence we find that P{r) is a second order quantity. This equation represents the force balance between pressure 
and gravity acting on the matter field. 

Taking it into account that P(r) is second order and expanding p as p = p^ 1 ' + p^ , the explicit junction conditions 
for the metric functions become 



dyB^ = dyC^ 



^P (J) , (at y = 0+). 



(2.23) 



The boundary condition obtained above is written in terms of the variables in the GN gauge. To interpret the 
conditions ( 2.23| ) in terms of the variables in the RS gauge, we consider the gauge transformation between these two 
gauges, which is defined by 



g ab (x) dx a dx b = g a b(x) dx a dx b , 



(2.24) 



with x a = x a + £, a (x). Since h yfJi and h yy vanish in both gauges, the infinitesimal gauge transformation £ a =£ a + £ a 
between them is restricted to the following form 



(J) 



(J) 



4a 2 



<i) 



(2.25) 
(2.26) 



U) (■;) 

Functions £ y (r) and £ r (r) are integration constants which appear as a result of y- integration. We discuss how these 
functions are determined soon later. 

We denote the difference between the metric in the RS gauge and that in the GN gauge as 



6A^(r,y) = A^(r,y)-A^(r,y). 



(2.27) 



The remaining metric functions 5B^ J ' and 8C^ are defined in the same way. Then, the gauge transformations are 
given by 



SA (J) = -- e Aty ? +A<p C 

t L ,w 



SB^ 

sew 



2 
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e +2 Cr 



B 



i ± i i ± i / ( 1 ) \ 2 ! : 

(i) ^ +0(1) r _( Cr ) +eW t (&) 



(2.28) 



2 < J > 2 < J > 
l r 



(J) 



I (1) \ 2 



Now we are ready to derive the equation that determines £ y (r). We eva luate the identity tj) tV — ip. y + Sip ty at y = 0. 
The expression for the left hand side is obtained by means of Eq.(2.17), while the right hand side is evaluated by 
substituting Eqs. (2.23), (2.25), (2.26) and ( [2.28 ). After tedious but straightforward computation, we obtain 



A P = 1lr {J) + e (J) 
* 6 



(2.29) 



with 



-1 t 1 ' 



S(r) = A -(e) 2 + / & B (1) (r',0)dr' - -(e) + f d r A £» +-(e r ) 2 + -(^ r ) 2 



(2.30) 
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The gauge freedom for the radial coordinate in the GN gauge has not been fixed. Although the simplest choice 



(J) 



might be to take £ r to vanish on the brane, for later convenience we impose the isotropic gaug e condition on the 
3-brane, B = C at y — 0. Rewriting this condition by using B and C with the substitution of ( 2.2S ), we obtain 



(i) 

r 



(2.31) 



where we have used (|2.8| ). Here we show the result only for the first order, because we do not need the explicit 



expression for £ r in the following discussion. 



(.J) 



Once we obtain an explicit expression for £ y , it is easy to d erive the boundary condition for the metric functions 
in the RS gauge from the j uncti on condition in the GN gauge (2.23). Especially, the boundary condition required to 
solve the master equation ( 2.19 ) is deduced by substituting the relations obtained above into the right hand side of 
the equation, A. y = A, y + SA, y . Imposing this boun dary condition is equivalent to add a 5-function source localized 
on the 3-brane. The 5-dimensional master equation (2.19) including the boundary condition becomes 



[C + U- X 8{y)] {a 2 A {J) ) = 2k 5 ^ J \r) S(y) + ha 2 ^ 1 ^^ - S 



(2.32) 



where 



SW(r) S « s -V |J, lrt 



P + K^dy(SA^)\y =0 



By using the Green function which satisfies 

[£ + U-H{y)] G(x, y; x', y') = S(y - y')5 3 (x - x') 
the formal solution for the master equation ( [2.32 ) is given by 

a 2 A {J) {r,y) = 2k 5 J dx' 3 G(x, y; x', 0) S (J) 

- 2e< J ' J d 3 x' J™ d y 'G( X ,y; X ',y 



(2.33) 



(2.34) 



S(r',y') + -„ /" a 2 (y")Q yy (r\y")dy 



(2.35) 



The factor 2 in the second term of ( 2.35| ) comes from Z 2 -symmetry. In §11] we discuss some basic aspects of the Green 
function, and solve the above equation at the linear order. In §[V we extend our analysis to the second order. 



III. GREEN FUNCTION AND FIRST ORDER PERTURBATIONS 

A. Green function 

We need the Green function to evaluate explicitly the formal solution ( 2.35| ). In the static case Green function is 

1 



G(y-,y,y-,y 



Att£R 



a 2 (y)a 2 (y')+£ / u m {y)u m {y')e' mH dm\ , 



(3.1) 



where u m {y) is the mode function, and we have introduced the notation 

R = |x-x'|. 

The explicit form of the mode function is given in terms of Bessel functions as 

u m (y) = N m [j 1 (me)Y 2 (m£/a) - Y x (ml)J % (ml t 'a)], 



(3.2) 



with N m = y/m£/ y / 2(Ji(m£) 2 + Yi(m£) 2 ). It is orthonormalized 



G 



dy 



2 / -rrUmUm' — S(m — w!) • 



(3.3) 



In particular, setting m! = 0, we have 



dy u m = (m 7^ 0). 



(3.4) 



The first term on the right hand side of (|]l]) is the contribution from the zero mode whose 4-dimensional mass 
eigenvalue is zero (to = 0). We denote this part of Green function by Go- The second term corresponds to the 
propagator due to the Kaluza-Klein states which have nonzero mass eigenvalues (m > 0). This term is denoted by 
G K . Thus, G = Gq + G k . 

For large separation R 3> i, the existence of the factor e~ rnR in the second term in ([Q|), i.e. Gk, implies that the 
integral over to is dominated by the contribution from small m. When the source is on the 3-brane (y' — 0), we can 
approximately evaluate Gk by expanding Bessel functions taking m£ as small, but here the Bessel function with the 
argument m£/a is to be kept unexpanded because a can be exponentially small. Then, the integration over m can be 
performed for the leading power of to explicitly. The result is given by H 



G(x,2/;x',0) 



a 3 {2a 2 R 2 + U 2 ) 
"^d{a 2 R 2 +£ 2 ) 3 / 2 



(3.5) 



For a small separation R £, the Green function is dominated by the contribution from modes with large to. In this 
limit, the Green function behaves as the ordinary 5-dimensional one w —1/Att 2 [R 2 + 1 2 {\ — a^ 1 ) 2 }. 
When we discuss second order perturbations later, we use the following inequalities 



0<-G(x,y;x',0) < 



a 2 (y) 

AttIR 2 



[R + l), 



(3.6) 



This inequality is suggested by the asymptotic form of the Green function, and it is confirmed by numerical calculations 
as shown in Fig{j]. 



log[-G] 




FIG. 1. We calculated numerically G and a 2 (R + t) / (A-kIR?) + G to show the inequalities of (3.6). Since there appears no 
logarithmic divergence, G and a 2 (R + 1) / (AixiR?) + G are positive definite at least in the region shown in the figures. 



B. linear perturbations 

We start with linear perturbations. From Eq.(^3E), is given by 

a 2 A (l) = ^ J dVG(x,y;x',0VW(r'). 



(3.7) 
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Suppose that is the radius of the star. When we consider the metric at a point outside the star r — £ or far from 
the 3-brane, we can safely replace G by the approximation (3.E) with the relative error of 0(£ 2 /r 2 ) ^. Furthermore, 
at a field point far from the star, the matter distribution can be repl aced with a point source. Then, we obtain 
A w {r,y) ~ ~4:G 4 Ma{2a 2 r 2 + 3£ 2 )/3(a 2 r 2 + £ 2 ) 3 / 2 . By using Eq. (fTlf) , we find fiW(r,y) - 4 G 4 Ma/3Va 2 r 2 + £ 2 . 

On the other hand, if we are interested in the met ric induced on the brane, we can also set y = 0. Even if we 
consider the interior of the star, the inequality ( |3.6| ) implies that the contribution to A^ from Gk is, at most, 
of 0((K 4 M*/r*)(£ 2 /r 2 )log(£/rA), where M* = 4tt / r 2 p{r)dr. Hence, it is small by a factor of 0{{i 2 /r 2 ) \og{£/rA) 
compared with the leading term. Then, neglecting the collections of this order or higher, we obtain 



B«(r,0) « -J-0 r A" 



3r 



where A 1 is the inverse of the Laplacian operator, and we have introduced Newtonian potential 

<t> = 47rG 4 A -1 /). 



These are transformed to the isotropic GN gauge by using ( 2.28]) with 



(i) 



t 1 ) 2 n 



which are derived from ( 2.29| ) and ( |2.3l| ), respectively. Finally, the resulting metric functions turn out to be 

- A {1) (r, 0) m (r, 0) = C (1) (r, 0) « -2#r), 
which agree with the result for the 4-dimensional Einstein gravity. 



(3.8) 
(3.9) 

(3.10) 

(3.11) 
(3.12) 



IV. SECOND ORDER PERTURBATIONS 



For later convenience, we quote the result obtained in the preceding section as 

2 f°° 

A (2) (r, 0) = A s + A s + A Q - - A" 1 J dy'a 2 Q yy , 



where 



A s 
As 

An 



2k 5 J dx /3 G(x,0;x',0)S (2) , 

' / dy'G(x,0;x',y') ^d r , {r' 2 A^(r',y')B^(r',y')) , 



2 / d 3 x 



4 / d 3 ~' 



■'" / dy' 

o 



G(x,0;x',y")cfc/" 



(4.1) 



(4.2) 



The source te rm I]( 2 ) for A^ is given in Eq.(2.33). The third term in (2.35), is separated into two pieces, Aq and the 
last term in ( |4.l| ), by performing integration by parts. To obtain the expression of this last term, we have also used 
J d 3 x' f^Gdy' = J d 3 x' J °° G dy' = (1/2)A -1 . Here the first equality follows from Eq. Q. 

Since we arc interested in the metric induced on the 3-brane, we also write down the expression for A™ at y = 0. 
Combining fl2~25| ), ( |2~28| ), ( f2~2S| ) and (p"!"]), we obtain 



A^ (r, 0) = A s + A s + A Q + ^ A-V 2 > + [ \ A' 1 H - \ (g.) - i« - A« f 



(4.3) 



s Here r* in the denominator is just inserted to adjust the dimensionality. It can be r instead of r*. 
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Here we note that the last term in (4.1) is canceled with a term which arises from the gauge transformation. 

The terms in the square brackets are evaluated just by substituting the estimate at the first order, and the first 
term As has structure similar to first order perturbations, i.e., the source term is localized on the brane. Thus, the 
evaluation of these terms is straightforward. What needs detailed investigation is the evaluation of As and Aq. Since 
the source terms of A$ and Aq, which are quadratic in the linear perturbation quantities, distribute through the 
5-dimensional bulk, it is necessary to evaluate a con voluti on of three 5-dimcnsional Green functions. 

Deferring the estimate of this convolution until § IV B , let us turn to the spatial components of the second order 

(2) 

metric perturbations. Although each spatial component depends on the choice of £ r , the gauge invariant combination 



- d r (rC^)\ = Ib™ - d r (rC^) - 8B { ^ + d r (rSC^)\ 

I ) y=0 I ) i 

= + ^d r SA^ - SB™ + d r (rSC^)\ 

I 2 ' 2 J y=0 



1 v=o 



dy Q 



(4.4) 



does not contain £ r . Here we have used Eqs. ( 2.10| ) and ( 2.1 8| ) with y = to eliminate B and C. Taking the isotropic 
gauge B^ = £( 2 \ the left hand side becomes — rB,r . Integrating this equation with respect to r, we obtain the 

(2) 

expression for B^(— C^). Here we note that the terms in the square brackets in Eq.(4.4) contains £ v , from which 
a term with y-integration, 

dy a 2 Q yy , 



arises. Combined with the last term with y-integration in Eq. (4.4), this term is reduced to the expression that does 
not contain y-integration. The detail of calculation is explained in Appendix A. Here we just quote the final result 



S^(r,0) = 4A^+A-^(^T^ + ^E+ 1 -Q yy -^(^) 2 )+Bs, (at y = 0) , 



(4.5) 



where 



By = - 



3 / ^ \ - 



(3^W +B«)S« - 1(3B« +A<$>) J r 2 ^dr 



y=o 



We find that the expression is reduced to a closed form written solely in terms of the quantities on the 3-brane except 
for^ 2 )(r,0). 



A. Recovery of the 4-dimensional Einstein gravity 

In this subsection, we evaluate the metric induced on the 3-brane at the leading order in f/r», and show that the 
result completely agrees with the one predicted by the 4-dimensional Einstein gravity. 
As for As, since the source term is localized on the 3-brane, we can approximate it as 

Ae«2k 4 A- 1 S^, (4.6) 

which is justified for the same reason explained in evaluating the induced metric at the linear order. 

As mentioned earlier, the point which needs careful analysis is the computation of ^4 5 and Aq. To evaluate these 
terms we need to evaluate a convolution of the 5-dimensional Green functions. Since G is composed of Go and Gk , the 
contribution from these terms is decomposed into several pieces depending on which combination of three propagators 
is used. For example, there is a mode coupling in which a zero mode propagator Go propagates the second order 
source produced by a product of KK mode contributions. We denote this mode coupling as 

[K,K;0], 

and similar labels are assigned for the other mode couplings, too. 

Here in the present subsection, we simply neglect the pieces containing the KK propagator Gk- Namely, we just 
take into account ^s[o,o;0] an( A -^q[o,0;0]- The justification of this approximation is given in the next subsections, 
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where we prove that the neglected pieces in As and Aq containing Gk are actually suppressed by a factor of 
0{{e/rl)\og{n/i)). 

As long as only the contributions from the zero mode are concerned, the first order metric functions are all constant 
in y. Since the source term of Aq contains differentiation of the first order metric functions with respect to y, we find 
that A, 



0[0,0;0] 



0. As for Agr 0i o;o]) after the y-integration, we find that it is reduced to 



(4.7) 



Now that t he eva luation of A*- 2 ' (r, 0) is straightforward. Substituting the evaluations of the first order quantities 
presented in § IIIB, we finally obtain 



and 



A™ (r, 0) « k 4 A- 1 (pW - 2( ^« + 3P) , 



£ (2) (r,0) = C^(r,0) « -A" 1 [re 4 (f> {2) - 2<£p (1) ) + (<£, r ) 



(4.8) 



(4.9) 



I t mu st be n oted t hat B x in (4.5) gives only higher order correction of 0(£ 2 k1) as it is easily shown by using ( 2.23 ), 
(2.28), (2.29) and (2.31) These results agree with those for the 4-dimensional Einstein gravity, whose brief derivation 
is given in Appendix B. 



B. Suppression of the KK-mode propagation 

In the preceding subsection the terms in As and Aq containing a KK-mode propagator were neglected. Here we 
show that the contribution from these terms is in fact negligible. We begin with discussing rat her general things. 
As for As, the terms As[o,o ; -K"]> As[o,i<";0] an d As[k,o-,o] vanish because of the orthogonality (3.4). Hence the terms 
to be investigated are 1.) A S [k,k-.o] > ^s[o,k-,k] ^-s\k,0;K] and As[k,k,k]- As for Aq, the situation is a little simpler. 
Recall the fact mentioned above that the zero mode contribution in the first order metric functions is y- independent. 
Thus we can say that Aq[ ^.,] and Aqu^;*] vanish. Therefore, all the terms that we need to consider for Aq are 2.) 
Aq[k,K;0] and Aq[k,K;K}- 



1. A 



S[K,K;0], As[K,0;K] and As[0,K;K] 



The source term of As is rewritten as 



1ft (^w) = l 



i r 



drr 2 A (1 



-ftA^ A (1 



Each source term has the form of 



fi(s) = / dm u m (0) e 
Jo 



The function /z(s) is bounded by 



< (Jt(s) < 



2s{s + £)' 



(4.10) 



(4.11) 



with Or — 1, ^rft, r 2 d 2 or r 3 / Q r dr r 2 . By using the orthonormal condition (3.3), we obtain the formula 
r°° _ i 

/ dj/a~ 2 G , A-(xi,y;x 2 ,0)Gx(x 3 ,y;x4,0) = — ^(fl 12 + R 34 ), (4.12) 

Jo 2(47r)^i?i 2 it 3 4 

where Rab = — Xs| and 



(4.13) 



(4.14) 
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We begin with As\k ,k-o] + ^-S[K,K;K 1 ■ We can derive an upper bound for the absolute value of each term as 



3 / dV I dy'G(x,Q;x',y') O^A^(n,y')O^A^(r 3 ,y') 



(i 



n=r 3 =r' 



™ / d 3 x' J°° dy'^l(R + l) 



8ir£r 2 



(4.15) 



- ri— r ; 



where .A^' is a part of A^ propagated by the KK mode propagator, and we have used ([0]). In below, we show that 
the each term on the right hand side of the above inequality is at most 0(£ 2 log(r*/^)). We express the respective 
term as 



1 r^ R +£ 



R? 



n=r 3 =r' 



with 



vir^r^^r 1 dya 2 A$(n)A£>(r 2 ) 
Jo 

K 2f pr, pr, pri+r-i /T 3 +r 4 

4 ' dr 2 r 2 i: (1 \r 2 ) / dr 4 r 4 £«(r 4 ) / dR 1 



(i 



2 / dR 34 n(R\2 + R34) ■ 

r3— r&\ 



(4.16) 



(4.17) 



where v(ri,r 3 ) has been rewritten by using Eq. ( p~12[ ). Introducing U{s) = - f* ds' ds"^(s"), 

k 2 i r* r* ( 

v(rx,r 3 ) = / dr 2 r 2 Y, {1 >{r 2 ) / dr 4 r 4 £W(r 4 ) h/(n + r 2 + r 3 + r 4 ) - U(\n - r 2 \ + r 3 + r 4 ) 

2^ir 3 Jo Jo ^ 

-Ufa + r 2 + |r 3 - r 4 |) + I7(|r a - r 2 | + |r 3 - r 4 |) ) . (4.18) 



By construction, U(0) = 0. The inequality (4.14) indicates that 



< U(s) < 



s log 1 



(4.19) 



As long as s £ r*, we find that J7(s) is at most 0(ilog(r*/^)). 

Let us consider the action of once or twice differentiation with respect to r% on U. Since r\ appears in the arguments 
of U only in the combination of r% + r 2 or |n — r 2 \, the differentiation with respect to r\ can be replaced with that to 
r 2 . Then, integrating by parts, finally the differentiation can be moved so that it acts on the source term r 2 T,^(r 2 ). 
As long as first or second derivative is concerned, differentiation of U does not appear as the boundary term. The 
same thing is true for the pair of r 3 and r 4 . Now, using the bound for U obtained above, we can conclude that 

0^}of}v{r\, r 3 ) is at most 0(£ 2 log(r„/^)) for ri,r 3 < r*. When r 4 ,r 3 > r», the argument of [i cannot be small. 

Then, we can use the bound [i{s) < £/2s 2 . Therefore, Q^}C$}v(n, r 3 ) is 0(£ 2 ) for r 4 ,r 3 > r*. 

To conclude that A S [k,k-,o] + As[k,k-,k] is at most 0((£ 2 /rl)\og(r */£)), we have to examine whether [0™ of} 
v i r ii r 3 )) ri=r3=r behaves well at r — * 00 and at r — * so that the operation of J d 3 x'R~ 2 (R + £) is well-defined. 

First we consider the large r limit. Since r 2 and r 4 are bounded by the radius of the star, [i{R\ 2 + i? 34 ) can be 
replaced with /x(ri + r 3 ) w £/{r\ + r 3 ) 2 for large r\ and/or r 3 , and hence v(ri,r 3 ) ss £ 2 {n 4 M*) 2 jr\r 3 {r\ + r 3 ) 2 . The 

operation of O^} and O^a does not make the fall-off worse. For small 7*1, 



v(n,r 3 ) 



K 2p pr, pr, r r 3+ r i 

41 1 dr 2 r 2 ^(r 2 ) / dr A r^^) 



r 3 Jo 
n 2 £ rr. 



r3-n\ 



dR 3 4 fi(r 2 + R34) + -jrp!\r 2 + R 3A ) + 0(r\) 



r 3 Jo 



dr? 



■±.,11/ 
6 

dr 4 r 4 £ (1) (r 4 ) / dR 34 fi(r 2 + R 34 ). 
J\r3— n\ 



(4.20) 



Therefore, 0^}v(ri, r 3 ) is finite at r 4 
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For As\k,0;K]i a similar expression is obtained as 

2 



with 



A-S[K,0;K] = ^ I dr ' r ' 



r\—r,r-2,—r^—r' 



(4.21) 



w(r 1 ,r 2 ,r 3 ) = r 2 I dft 2 I dy' G K (*i, 0; x 2 , y^A^ (r 3 , y') 



« /-r, rr 1 +r 2 rr 3 +r4 

— =- / dr 4 r 4 E (1 >(r 4 ) / di? 12 / dR 3A fx(R 12 + R 3 a) , 

4rir 3 Jq J\ri-r 2 \ J\r 3 -n\ 



(4.22) 



and ^Ig 1 - 1 represents the part of .A^ 1 ) propagated by Go- In the same way, we can show that -kt] is at most 

0((* 2 /r?)log(r,/*)). 



2- -Aq[K,K;0] an ^ ^-Q[K,K;K] 



Aq[k,k-,o] and Aq^^.^j are also bounded as 

|^Q[K,K;0] + ^•Q[iir,Jr; J K]| 









R + l 


< 




f d 3 x' 


R 2 


2-kIJ 








R + £ 


< 




f d 3 x' 




8tt£ J 


R 2 



dya i (y)Q vv (x, y) 



y=o 



dy ((aUW),^ + (a 4 B«), tf Bg) + 2(a 4 C«),,4 1) ) j. (4.23) 

The first term in the curly brackets is at most 0{{£ 3 /rf) log(r*/-£)). By using the relations obtained from Eq. ( 2.32 ), 

-MA^.y = -a- 2 AA^ + 2 K5 £ (1) <%), (4.24) 
a 

i(a 4 S«),, = a~ 2 -d r A^ - ^ [ drr 2 £«%), (4.25) 
a 4 r r J J 



and = — (^4^ 1 - ) + B^)/2, the second term is reduced to the form in (4.16). Therefore, the contribution to Aq 
from this term is at most 0((£ 2 /r 2 ) log(r*/£)) for the same reason. 



V. SUMMARY 

In this paper we developed second order perturbations in the RS single brane model, restricting the configuration 
to the static axisymmetric one. From the 5-dimensional Einstein equations, we derived the master equations for 
second order perturbations. At the level of linear perturbations, we can use the RS gauge, in which all perturbation 
equations are decoupled. Since the transverse-traceless condition cannot be imposed on the second order, the second 
order perturbative equations are inevitably coupled. As we have shown, however, the 4-dimensional spatial components 
of the second order metric function B^ and turned out to be concisely represented by the temporal component 
A^ 2 ' with the first order metric functions. Therefore, the problem was reduced to solving a single differential equation 
for A^ . Once we solve for A^ 2 \ the other metric functions B^ and C^ 1 follow from it. Further, to discuss the metric 
induced on the 3-brane, we introduced the GN gauge in which a hypersurface with constant 5th-coordinate coincides 
with the location of the 3-brane. we gave the second order gauge transformations between the RS gauge and the GN 
gauge explicitly. 

Based on this formulation, we first discussed the zero mode truncation for second order perturbations. It was shown 
that the metric induced on the 3-brane evaluated by using the approximation of the zero mode truncation exactly 
agrees with that for the 4-dimensional Einstein gravity. 

Next, we evaluated the contribution to the metric functions from the KK modes. Since the mode-by-mode analysis 
shows a pathological feature even at the level of linear perturbations, it is necessary to sum up all the mass eigenvalues 
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to handle the KK mode interactions. Performing such an analysis, we have confirmed that the correction due to the 
KK mode coupling on the induced metric is suppressed by a factor of 0{{£ 2 /r 2 ) log(r*/£)), and there appears no 
pathological behavior. We therefore conclude that second order perturbations in the RS single brane model behave 
well and the result basically agrees with the prediction by the 4-dimensional Einstein gravity. The relative order 
of the correction is 0((£ 2 /r 2 ) log(r»/^)). In the language of the post-Newtonian (PN) analysis of the 4-dimensional 
Einstein gravity, the order of A^ and B^ is 1PN, and that of B^ is 2PN. According to the parameterized post- 
Newtonian (PPN) formalism, the PPN parameters at the 1PN order, (3 and 7, are observationally constrained to the 
accuracy of about 0.1% or so |p2| . Although this accuracy might be improved in the future, there will be no chance, 
unfortunately, to find the evidence for the large extra-dimension from the precision test of the solar system since the 
predicted deviation in (3, 7 is 0((£ 2 /r 2 ) \og{r*/£)). 

Of course, since our discussion developed in the present paper is restricted to second order perturbations with 
spherical symmetry, the gravity induced on the 3-brane in the RS single brane model might deviate from the 4- 
dimcnsional Einstein gravity in more general situation. Hence, it would be necessary to develop a more complete 
proof of the coincidence including all higher order terms would. Although we did not discuss the RS two brane model, 
it is also interesting to study the post Newtonian correction in this model. The analysis is now in progress and will 
be reported soon in a separate paper. 
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APPENDIX A: DERIVATION OF S (2 



As was mentioned in the text below Eq. (T4), the expression given in the last line of Eq. (^J) contain the terms 
with y-integration. They are explicitly written as 



A- 1 



dya 2 Q 



uii 2 &r I dy T Qyr 



(Al) 



The first integral comes from £ y contained in SA^ 2 \ SB^ and SC^. In obtaining the expression (4.5) from 

how to reduce these terms will be the only non-trivial manipulation. The rest of computation is slightly complicated 

but almost straightforward. 

To rewrite the first integral, we use the relation 



~Qyy(r t 0) = \ jf° dy dy {a^A^f + {a'B^f + 2(a^) 2 }} 



3 ,00 

tJ 



1 f°° 1 

dya 2 Q yy + 5 / dy-[(3A^+B^)(aU^), y + (3S« + A^)(a*B^), v ], 

o Jo a 



where we have used the traceless condition = — (^l^ 1 ^ + B^)/2. On the other hand, the second part is rewritten 
as 



1 



1 



-ird r I dyr Q yr = —-Q 



dy [(3^W + B$)AAW + (30W + A^)AB^}, 



where again we have used the traceless condition. After the substitution of these relations and ( 4.24 ), the remaining 
y-integration is just 



A 



We can derive 



from Eq. (4.25), and (3S^ + A$) = —rB)ry from Eq. (2. IS). With the aid of these relations, it is easy to see that 
the remaining y-integration (|A2|) can be performed. 



-<9 y a 4 <9 y 



AB« + -B«)-^ / drr 2 ^5{y), 



(A2) 



(A3) 



2(1) 
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APPENDIX B: SECOND ORDER PERTURBATIONS IN THE 4D EINSTEIN GRAVITY 



In this Appendix, we give second order perturbations in the 4-dimensional Einstein gravity for the comparison with 
the zero mode truncation of the gravity in the RS single brane model. Although the results are well known as the 2nd 
post-Newtonian analysis |33| , we present a brief derivation for the following two reasons. 1) Since we are working in 
a specific gauge, we need to consider a gauge transformation to compare our results with the expression presented in 
a different gauge. 2) To compute metric perturbations in our restricted situation from the beginning is much easier 
than to follow a literature in which unrestricted cases are discussed. 

We assume that the 4-dimensional metric is static and isotropic, 



ds 2 



Up to the second order, the 4-dimensional Einstein equations with the energy momentum tensor given in (2.2) are 



Ri = - 2 (B 



l)AA--A r (A + B) 



«4 



rt' r 



(B-l) (A(A + 2B)--(A + B) 



(P + 3P), 

+ -A_ r (B - A) r 



fi'4 



(p-P). 



(B- 1) AB + -(A + B) 



-BJA + B) 



K 4 



(P-P)- 



(Bl) 
(B2) 
(B3) 



We expand A and B to the second order as fl2.9| ) . So lving these equations at the linear order, we obtain A^ — 
— P/ 1 ) — 2<j), where cj>(r) is similarly defined as in ( |3.10 ). Putting these results into Eq. flBlf ), the equation for A 1 ^ 1 
becomes 



AA {2) = k 4 (/0 (2) + 3P) - 40A</> , 
which agrees with flQj|). Eliminating AA^ and P from Eqs. pi] ), Q and (|B3|), we obtain 

AB< 2 ) 



-K 4 p 



(2) 



bAc 



(B4) 



(B5) 



Hence it was shown that the result obtained by the zero mode truncation ([O]) agrees with that for the 4-dimensional 
Einstein gravity (B5). 
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